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The technique of truncated moments of parton distributions allows us to study scaling 

Q I violations without making any assumption on the shape of parton distributions. The numerical 

implementation of the method is however difficult, since the evolution equations for truncated 

moments are not diagonal. We present a simple way to improve the efficiency of the numerical 

solution of the evolution equations for truncated moments. As a result, the number of truncated 

, moments needed to achieve the required precision in the evolution is significantly smaller than 

qj I in the original formulation of the technique. The method presented here can also be used to 

i-^ ' obtain the value of parton distributions in terms of truncated moments, and therefore it can 

be viewed as a technique for the solution of the Altarelli-Parisi equations. 
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The measurement of deep-inelastic scattering structure functions is a central issue in strong 
interaction physics: it allows the extraction of the parton distributions of hadrons, which, 
though in principle computable, are determined by the nonperturbative dynamics of the theory, 
and must be treated as unknown phenomenological parameters. A detailed understanding of 
these quantities is an essential ingredient of phenomenology at hadron colliders [|l], 0] . Further- 
more, the measurement of scaling violations allows a determination of the only free parameters 
in the QCD lagrangian, the strong coupling a^, and of the gluon density in the nucleon. 

As is well known, scaling violations of parton distribution functions are described by the 
Altarelli-Parisi (AP) evolution equations 0. These are integro-differential equations whose ker- 
nels are presently known to next-to-leading order [Q (for a handful of operators, the anomalous 
dimensions are known to next-to- next-to-leading accuracy 0). The coefficient functions that 
determine the relation between measured structure functions and parton distributions have also 
been computed up to next-to- next-to-leading order 0. 

There are different techniques to solve the Altarelli-Parisi equations, and the choice among 
them depends on the kind of problem at hand. Usually we solve the evolution equations 
analytically by taking their Mellin transform, which turns convolution products into ordinary 
ones, and therefore the x-space integro-differential equation into a set of independent ordinary 
first order differential equations. A parametrization of the distributions is assigned at some 
initial scale, and the parameters are then determined by fitting to data the evolved distributions. 
Mellin moments of structure functions, however, cannot be measured even indirectly, since they 
are defined as integrals over the whole range < a; < 1, and thus require knowledge of the 
structure functions for arbitrarily small x, i.e. arbitrarily large energy. 

We can solve this problem using the Altarelli-Parisi equation to evolve parton distributions 
directly: the scale dependence of any parton distribution at xo is then determined by knowledge 
of parton distributions for all x > Xq, i.e., parton evolution is causal. In fact, through a 
judicious choice of factorization scheme 0, |^ all parton distributions can be identified with 
physical observables, and it is then possible to use the Altarelli-Parisi equations to express the 
scaling violations of structure functions entirely in terms of physically observable quantities. It 
is, however, hard to measure local scaling violations of structure functions in all the relevant 
processes: in practice, a detailed comparison with the data requires the solution of the evolution 
equations. 

As pointed out, the solution of the evolution equations requires an assumption on the x 
dependence of the parton distributions at the initial scale; a frequently-adopted input is for 
example p 

q{x, Ql) = ao x'^Ml - ^r P{^: , as, . . .) , (1) 

where Qq is a reference scale. The parameter ai is associated with the small-x behavior while 
02 is associated with the large-x valence counting rules. The term P(x; , as, . . .) is a suitably 
chosen smooth function, depending on one or more parameters, that adds more fiexibility 
to the parton distributions parametrization. It has however become increasingly clear that 
in practice this procedure introduces a potentially large theoretical bias, whose size is very 



hard to assess [0]. In ref. |Ty] it was proposed to adopt a functional method to keep this 
theoretical error under control. Another suitable way to minimize the bias introduced by the 
parton distributions parametrization is to project the parton distributions on an optimized 
basis of orthogonal functions. Different methods have been suggested with suitable families of 



orthogonal polynomials {e.g. Bernstein [11], Jacobi [12| or Laguerre polynomials [^]) as basis 
of function. 



A different approach has been suggested in refs. [TJ, 0, which makes use of truncated 
moments of parton distributions. Truncated moments are defined in analogy with ordinary 
moments, but the integration in the Bjorken variable x is now restricted to the subset xq < x < 1 
of the allowed kinematic range < x < 1. As a consequence, the corresponding evolution 
equations are not in diagonal form: the evolution of the moment of order no depends on all 
moments of order UQ+k, with A; > 0. The solution of the evolution equations is more complicated 
than in the case of ordinary moments, but it can be performed by taking only a finite set of 
truncated moments. Assigning their values at a reference scale as input parameters, we can fit 
the value of as or the first moment of the gluon; the initial values of truncated moments can 
be obtained directly by data. In this way, no assumptions are made on the shape of parton 
distributions. The fit is only affected by the experimental errors and the theoretical uncertainty 



that affects the Q^ evolution can be easily kept under control. However, it was shown in ref. [15 



that the number M of truncated moments needed to achieve a precision on the evolution of the 
lowest moment comparable to that of other techniques is rather large (M ~ 150), and in some 
cases it may lead to problems in the numerical implementation of the method. This problem 
was overcome in ref. |1^ by showing that, in practice, it is sufficient to parametrize the parton 



distributions using the first few (between 7 and 10) truncated moments, plus the value of the 
parton distributions at x = a;o. 

In this paper we present a different way to improve the numerical efficiency of the method of 
truncated moments. We will find that integrating by parts the RHS of the evolution equations 
for truncated moments allows us to use a significantly smaller number of truncated moments 
involved in the evolution (M ~ 10). We will see that the integration by parts introduces an 
explicit dependence on q{xo, Q"^) of the AP equations, which complicates the solution of the 
evolution equations for truncated moments. We will show how to circumvent this difficulty. 
As a by-product, a formula is derived for the evolved parton distribution at all values of x 
larger than Xq, in terms of the first M truncated moments; the method provides therefore an 
alternative way of solving the evolution equations. 

We begin by studying the unpolarized non-singlet case at leading order, leaving at the end 
the extension to next-to- leading order. The Q^ dependence of parton distributions q{x, Q'^) is 
governed by the Altarelli-Parisi equations [§] 

^,,M'-) = ^r-^p('^-,a.(Q-,),(y,Q'), (2) 

at zir Jx y \y J 

where t = logQ^/A"^. The evolution kernels P(a;, ^^(Q^)) are perturbatively computable as 
power series in as- In the non-singlet case, q{x,Q'^) is simply one of the fiavor non-singlet 
combinations of quark distributions and P(x, as(Q^)) the corresponding splitting function. 
The truncated moments of a generic function f{x) are defined as 

fn{x,)= f'dxX^-'f{x). (3) 

Jxo 



The evolution equations for truncated moments of parton distributions were derived in [|14|, |15 . 



One finds that the truncated moments of q{x, Q"^) obey the equation 

J X 

asm 



with 
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dt' a{t') ; a{t) 
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In ||14| , |15[ it was shown that eq. (^) can be written as 



d A/+no 

-l-qniXo,Q'^) = Yl CniqiiXo,Q'^) 
l=no 



(4) 

(5) 
(6) 

(7) 



where Cni are the elements of a triangular matrix, and only a finite number M of truncated 
moments is taken into account. It was also shown that in order to reach a precision of 5% on the 
RHS of the evolution equation for the lowest value of no a large value of M is needed: M ~ 150. 
This makes it difficult to solve the evolution equation as we need a large numerical precision. 
In the following we will show how these difficulties can be overcome. For later convenience we 
will set no = 1. 

We now integrate by parts the RHS of eq. (§). We get 



xo 



dyy--\iy,Q')a 
^0 \y 



d 



Gn{xo,y)y''-'q{y,Q')] - dyGn{xo,y)- (y''-'q{y,Q')) 



where 



Gn{xQ,y) = dzGni — 

J xo \ Z 



(9) 



(the lower integration bound is irrelevant here; it has been chosen equal to Xq for later conve- 
nience). Using the definition of Gn{xQ,y) and eq. (j^), we get 



Gn{xo,y) = dz dxx"' ^P{x 

J xo Jxo/z 



dxx^'-^Pix) ^^ '- -^ ^^° 



dz = yGr. 



xo/y Jxo/x \ y 

By taking the Taylor expansion of Gn{xQ,y) around y = 1, we obtain 



— XqCj 



Xo 
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yqn{xo,Q^)= Gn{xo,y)y'' ^q{y,Q^) -Y. 

CLT '- -'^0 n p'. Jxo 



p=0 



_ dy{y-lY^(y^~'q{y,Q' 

^0 :-^ v\ Jxo dy ^ 



(10) 
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where 



9i{.xo) 



dP ~ 

- — Gnixo, y) 
dyP ^ ^ 



y=l 






9'n-\^0] 



(12) 



y=l 



The functions Gn{xo, y) are regular in the whole interval [xq, 1]. In fact, the Gn{xo/y) are regular 
for all values of y except y = Xq, as they contain singular terms proportional to log(l — Xo/y). 
However, these terms are integrable, and independent of n. Thus, Gn{xo,y) is regular in the 
limit y ^ xq and tends to zero. Furthermore, we observe that the Taylor coefficient of order p 
of Gn{xo, y) is equal to that of Gn{xo/y), times a factor 1/p (see eq. (|12D ). For this reason, the 
convergence of the expansion of Gn{xo,y) is faster than that of Gn{xo/y). 
Integrating by parts the second term of the RHS of eq. ([TT|), we have: 



d_ 



qn{xo,Q^ 



Gnixo,y)y''~My,Q') 



xo 
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xo 



q{y,Q' 



(13) 



Truncating the series, expanding the binomial [y — ly ^ and imposing g(l, Q^) = (this is our 
only assumption on the behavior of the parton distributions), we get 



dr 



M 



qn{,XQ,Q^ 



X, 



n-1 



q{xo,Q^) Y. 



p=Q 



fc=0 



p\ 



xo)qn+k{xQ,Q^ 



(14) 



where 



M 



Cnk (Xo) = Y. 



p=k 



-l)P+%"(Xo) 
k\{p-k)\ 



(15) 



Defining the triangular matrix 




ciy.-r^ ii>k) 





(/< A;) 
we can finally write the truncated evolution equation as 



(16) 



-i-qn{xo,Q'^) 



M 



„n-l 



q{xo,Q^) Y 



9?,ixo) 



M 



p=0 



Xo - I)'' + V Cniqi{xo, Q^ 



;i7) 



Notice that the first term in the RHS of eq. (p^ vanishes in the limit M — > oo and the original 
expression given in refs. O, |T3| is recovered. However, for finite values of M this term must be 
taken into account (in a sense, it is the price we have to pay for the better convergence of the 
expansion after the integration by parts). This term poses special problems because it depends 
on the value of the parton distributions at x = xq. In the following we will show how to obtain 
an approximated expression of q{xo,Q'^) in terms of a finite number A^ (not necessarily equal 
to M) of truncated moments. The evolution equation ( [T7| ) will then be solved with a technique 
similar to that presented in [jl4|, p!5| . 

We begin by taking the Taylor expansion of g(x, Q"^) around x = yo: 



qix,Q') = J2VkiQ')i 



X 



yo 



,fc-i 



fc=i 



The initial point of the expansion, yo, must be carefully chosen. Parton distributions parametri- 
zed as in eq. (|T]) are non-analytical in a; = 1 when the exponent 02 is not an integer; and even in 
that case, an essential singularity in a; = 1 is generated by perturbative evolution. One should 
therefore choose yo < (1 + xo)/2, so that the expansion ( ]T8| ) is convergent everywhere in [xq, 1). 
The series will not be convergent in x = 1, no matter what yo is; however, the singularity in 
X = 1 is integrable, and the term-by-term integration is allowed using the Lebesgue definition 



of the integral (see ref. [|T5| for details). We have therefore 

•^^■0 k=i 



where 



Pjk{xo,yo) = dxx^ ^ 

Jxo 



X - yo) 



\k-l 



(19) 



(20) 



Our task is now to find a way of inverting eq. (|T9D, in order to express the coefficients ^fe(Q^) 
in terms of the truncated moments qj{xo,Q'^). This can be done in the following way. Define 
a matrix (3^^ by 



Phi 



otherwise 



(21) 



where /3*^^-* is the N x N upper left submatrix of /5. For example, in the case N = 2 the matrix 
P~^ is such that 



p-'-P 



/I n /^13/322— /323/3l2 

' ^ ^ det/3(2) 

fl 1 fea/^ii-Aafei 

^ ^ det/3(2) 





...\ 



•/ 



(22) 



Multiplying eq. ([T9|) by /? ^ on the right, we obtain 



T.Pij' ^A^o,Q') = T.T.Pi,' P^kVkiQ') 



(23) 



i=i 



k=ij=i 



where for simplicity we have not shown the dependence on xo,yo- Using the definition oi (3 ^ 
we get 



JV _ 00 AT _ 

j=l k=N+lj=l 



(24) 



for i < N. Substituting eq. (p^ ) in eq. (jT^) gives 



N N 



q{xo, Q^) = E E/^fc/ ?i(^o, Q^) {xo - yo)'' ^ + R{xo, yo, Q^ 



(25) 



k=ij=i 



where 



Rixo, yo, Q"' 



k=N+l 



[xo - yo 



,fc-l 



N 

E(^o 



TV 



yo) 



"EA^^A 



'jk 



i=i 



(26) 



We have thus obtained an approximate expression of q{xo, Q"^) as a function of the first N 
truncated moments of g, eq. p5|); the quantity R in eq. (pB]) represents the error on this 
reconstruction. The quantity in square brackets in eq. (|26|) is independent of the parton dis- 
tributions, and can be computed for any A^ and k starting from the coefficients j3ij, given by 
eq. (pO|) . The analytic expression of this quantity is very comphcated. We have checked that, 
for yo = (1 + a^o)/2, it decreases as [(xq — l)/2]'^~^, for any value of A^. Therefore, -R(xo, yo, Q^) 
vanishes, for A^ — *■ oo, at least as fast as the remainder of order A^ of the Taylor expansion in 
eq. (0). 

In order to assess the accuracy of our approximation, we have computed the percentage 
error given by ratio |i?/g(xo, Q'^)\ for some representative choices of the parton density, namely 



g(a;,Q2 



X 



112 



with 02 = 2.5,3.5,4.5. We have fixed xq = 0.1 and yo = {1 + xo)/2. 



The results are shown in Table |l|. We see that an excellent approximation is achieved already 
with N = 5, independently of the value of the large-x exponent 02. The accuracy increases 
with increasing A^; however, it should be noted that a numerical evaluation of the matrix j3~^ 
requires a numerical precision which also rapidly increases with A^. Therefore, for a practical 
implementation of the method, A^ cannot be very large. We see from Table |l| that for 5 < A^ < 
10 the accuracy is already better than 10~^ in the cases we have studied. We conclude that 



xo = 0.1 


N 


02 = 2.5 


02 = 3.5 


02 = 4.5 


5 


3.3 X 10^^ 


3.2 X 10"^ 


1.4 X 10-3 


10 


3.8 X 10-6 


5.4 X 10"^ 


1.4 X 10-^ 


15 


3.2 X 10-^ 


1.8 X 10-s 


1.8 X 10-9 



Table 1: Precision in the reconstruction of q{xo,Q^) = (1 — Xq)""^ in terms of a finite number 
N of truncated moments, for different values of N and for three different choices 0/02. 



N [■ N 



q{xo,Q^) ~ E 



i=i Lfc=i 



J2l^klixo,yo)ixo-yo 



,fc-i 



qj{xo,Q'^) 



d ^^ ^^ 

-i-qn{xQ, Q^) = E Cni qi{xo, Q^) + E ^'ni^ 11(^0, Q^ 
"•"^ 1=1 1=1 



(27) 



to an accuracy of about 10-^ for A^ = 5, independent of the detailed shape of q{x,Q^), and 
rapidly increasing with A^. 

We are now ready to re- write the original evolution equation ([T7| ) using our result eq. { \I7\) . 
We get 



(2^ 



xo = 0.1 


M 


niM 


'^l,M,6 


"^2,1/ 


'^2,A/,6 


5 


0.62 


0.14 


0.14 


0.020 


10 


0.48 


0.07 


0.12 


0.016 


20 


0.33 


0.03 


0.09 


0.009 


40 


0.20 


0.01 


0.05 


0.004 



Table 2: Comparison between percentage errors for the first and the second truncated moment 
at LO with N = Q, Xo = 0.1, y^ = {I + a;o)/2 and q{x, Q^) = (1 - xf-^. 



where Cni is defined in eq. 



and 






n-l 



M 



i:M^(x„ 



p=0 



pi 



N 



Lk=l 



yk-1 



(29) 



We now turn to a test of tlie accuracy of the evolution equation. We will also compare 
the accuracy achieved with the method presented here, and that of refs. |T^, |1^. The original 
evolution equation (|) and its truncated version, eq. (j^), can be written schematically as 



A. 

dr 



qn{Xo,Q'^) = Sn 



d_ 
d^ 



gn(Xo,Q^)=^r 



(30) 



respectively. Therefore, the quantity 






Q(M) 



(31) 



is a measure of the accuracy of the method adopted in refs. 0, 0. Similarly, we write eqs. (^Sf ) 
in the form 



^gn(xo,Q^) = 5r-^)+Tr'^) 



and define 



n 



b 

n,M,N 



1- 



Q{M~1) I rp{M,N) 
n '^ n 



(32) 



(33) 



a,b 



to test the error of the method presented above. The values of 7?,„'^, computed at leading order 
with xo = 0.1 for n = 1 and n = 2, are shown in Table ^ for different values of M and A^ = 6. 
We observe that the error TZ'^ j^j ^ computed with the technique presented here is always much 
smaller than the corresponding error of refs. [|l^, TZnM- ^^ accuracy of less than 10% can be 
achieved with a relatively small value of M. 



The complete solution of the evolution equation, LO and NLO terms, can be written as [T^ 



qn{xo,Q^) = R ^ 



ol'^ 



+ ae 



(7+60 )-r 



rfae-(^+^«)"(($i + I)i)e^° 



Rqn{xo,Ql) 



(34) 



with the initial condition 

Qnixo, Ql) = q'^\xo,Ql) + a{0)ql^\xo,Ql) , (35) 

and 

R{Co + Do)R-' = diag(7i,...,7Af) = 7, 
Ci + Di = R{Ci + Di)R-\ 

The matrix R that diagonahzes Cq + -Do must be computed numerically. This is not a problem, 
since, as we have seen, its dimension M does not need to be too large. 

In conclusion, we have shown that the evolution of truncated moments can be computed, 
to a degree of accuracy sufficient for practical purposes, by solving a system of a reasonably 
small (~ 10) number of coupled differential equations. The improvement of the numerical 
efficiency presented here can be straightforwardly extended to the unpolarized singlet case, and 
to polarized partons as well. The tests we presented are only for the LO equations. We have 
checked that the inclusion of NLO terms does not modify our conclusions. The technique of 
truncated moments can now be easily implemented numerically for phenomenological purposes. 

We have also shown that truncated moments provide, through eq. (|27|), a parametrization 
of the parton distribution itself, which has the advantage of being free of theoretical biases on 
the shape of the distribution at a given scale. 
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